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Abstract We analyse whether hierarchical formation models based on A cold dark 
matter cosmology can produce enough massive red galaxies to match observations. 
For this purpose, we compare with observations the predictions from two pub- 
lished models for the abundance and redshift distribution of Extremely Red Objects 
(EROs), which are red, massive galaxies observed at z > 1 . One of the models in- 
vokes a "superwind" to regulate star formation in massive haloes and the other sup- 
presses cooling through "radio-mode" AGN feedback. The first one underestimates 
the number counts of EROs by an order of magnitude, whereas the radio-mode AGN 
feedback model gives excellent agreement with the number counts of EROs and red- 
shift distribution of K-selected galaxies. This study highlights the need to consider 
AGN feedback in order to understand the formation and evolution of massive galax- 
ies at z > 1 . 



1 Introduction 

In hierarchical models, dark matter structures grow hierarchically, i.e., more massive 
haloes form most recently. However, this does not necessarily imply that the most 
massive galaxies formed the most recently. If the formation of galaxy mass tracked 
that of the dark halo, then there would be a problem for the formation of massive 
galaxies seen at z > 1 . Extremely Red Objects (EROs) allows us to study this point 
since they are red - (R — K) > 5 , massive galaxies, observed at < z >~ 1.1 L6| (this 
value depends on the magnitude limit of the observational survey). The extremely 
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red optical-near infrared colours of EROs, indicate either a starforming galaxy with 
heavy obscuration or an aged stellar population with little recent star formation |fT9ll . 
In the latter case, the implied age of the stellar population is uncomfortably close to 
the age of the Universe at the redshift of observation |20|. 

Reproducing the abundance of EROs has previously eluded many hierarchical 
models. Recently there has been much development in the modelling of the forma- 
tion of massive galaxies. Typically, hierarchical models produce too many massive 
galaxies at the present day unless some physical process is invoked to restrict gas 
cooling in massive haloes (see e.g. [21 for a review of plausible mechanisms). The 
most promising candidates are the heating of gas in massive haloes by "AGN feed- 
back" Q and the ejection of gas from intermediate mass haloes in a wind, which 
could be driven by supernovae 1 1 1 or by a QSO phase of accretion onto a black hole 

m. 

The bands used here are centred at 0.65 flm, R band, and 2.2 fim, K band. All 
magnitudes are on the Vega system, unless otherwise specified. 

The results described in this proceedings are taken from Gonzalez-Perez (2008) 

m. 



2 Galaxy formation models 

We study the abundance and properties of EROs in a ACDM universe using the 
semi-analytical galaxy formation code GALFORM 17|. Semi-analytical models ap- 
ply simple physical recipes to determine the fate of the baryonic component of the 
universe. 

For this work we focus on two published developments of GALFORM fl]: the 
Baugh et al. 1 1 1 and Bower et al. 1 3 1 models. We extract predictions for EROs galax- 
ies without adjusting the values of any of the model parameters. These parameters 
were fixed with reference to a subset of observational data, mostly at low redshift. 
Both models have some notable successes at higher redshifts: the Baugh et al.O] 
model reproduces the number counts of sub-mm selected galaxies and the luminos- 
ity function of Lyman-break galaxies, whereas Bower et al.|3| model matches the 
evolution of the K-band luminosity function and the inferred evolution of the stellar 
mass function. 

Table[T]lists the cosmological parameters used in each model and the main differ- 
ences between them. For our study, the key difference is the physics used to suppress 
star formation within bright galaxies: Baugh et al.|[T| introduce an outflow wind that 
ejects cooled gas from haloes, while in Bower et al.|f3l| model AGN feedback pre- 
vents gas from cooling in massive enough galaxies. 
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Table 1 Basic differences between the two models used in this work: 



Baugh et al. 2005 Bower et al. 2006 





0.3 


0.25 


Ao 


0.7 


0.75 


Qb 


0.04 


0.045 


On 


0.93 


0.90 


h 


0.7 


0.73 


Dark matter halo merger trees 


Monte Carlo 


N-body 


Quenching of star formation in bright galaxies 


Superwind 


AGN feedback 


Time scale for quiescent star fomiation 


Independent of time 


Dependent on time 


Bursts triggered by 


Mergers 


Mergers and disk instabilities 


Burst Initial Mass Function 


Top heavy 


Solar nei ghbourhoodi, 1 4 1 


3 Results 



In Fig. [U we compare the observed EROs number counts with the predicted ones, 
for different colour cuts defining the ERO population. The predictions of the Bower 
et al.|3 1 model reproduce the observed abundance of EROs impressively well for all 
the colour cuts shown, which is remarkable as none of the model parameters have 
been tuned to achieve this level of agreement. 

The agreement is less impressive for the Baugh et al.fT| model. Fig.[T]shows that 
this model typically underestimates the number counts of EROs by more than an 
order of magnitude. This lack of EROs cannot be blamed on a fluctuation in the 
observed counts due to sampling variations arising from the small fields and strong 
clustering of EROs. 

Some fraction of EROs are likely to be heavily dust extincted. It is important 
therefore to make a realistic calculation of the degree to which galaxy colours are 
extincted. We have experimented with changing and improving the dust treatment, 
without finding any significant change in the predicted ERO counts. We also found 
out that shortening the burst duration in Baugh et al.ll] model had little impact 
on the predicted number of EROs. The small number of dusty EROs predicted by 
this model is remarkable, and this could in part explain the under-prediction of the 
total ERO counts. Therefore, the disagreement between Baugh et al.| 1 1 model and 
observations suggests that the superwind feedback used in this model could be too 
efficient for bright galaxies at z > 1, delaying their formation. 

The predicted differential redshift distribution of galaxies brighter than K — 2Q 
and EROs selected with different colour cuts are shown in the left panel of Fig. ID 
It can be seen that Bower et al.|3| makes an accurate prediction of the redshift dis- 
tribution of K-selected galaxies up to z ^ 1.3. Such a good agreement is in contrast 
with the Baugh et al.||T| model which predicts a tighter redshift distribution around 
a lower redshift value. 

Fig. |2] shows that the observations of significant numbers of EROs with z > 1.5 
favours the Bower et al. 1 3 1 model; in the Baugh et al. [IJ model it is much less likely 
that one would find galaxies at such redshifts. 
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Within the UKIDSS survey data IfTSi a clear tendency was seen for redder EROs 
to be at higher redshifts. The same tendency is found for both models, as seen in the 
left panel of Fig. |2] 

For further comparison, we plot in the right panel of Fig. |2]the observed me- 
dian redshift distribution for both the K and R bands and compare with the models 
predictions. The Bower et al.|3| prediction is consistent with the observations in 
both bands. The Baugh et al.[lj model matches the data better in the R-band than it 
does in the K-band. This model underestimates the median redshift for K-selected 
samples faintwards of K 16. 




Fig. 1 K-band cumulative number counts for EROs selected by their {R — K) colour: {R — K) > 5 
in top left panel, {R — K) > 5.3 in top right panel, (S — A") > 6 in bottom left panel and {R — 
A") > 7 in the bottom right panel. The solid lines corresponds to the predictions from Bower et 
al.O and the dashed ones to the Baugh et al.[T] model. The source of the observational data 
is given in the legend fT0lT7)'20)fT9)'22l|l[T8l[l5l[l6l of each graph and it can be find at: 
[http://segre.ieec.uab.es/violeta/rawEROSdata .txt_ The errors shown are Poisson. 
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Fig. 2 Left The redshift distribution of EROs with K < 20, defined by {R — K) colour redder than 
5 (green histogram) and 7 (pink). The black histogram represents the redshift distribution for all 
galaxies with K < 20. The marks on the x-axis show the median of the distribution plotted with the 
same colour. Observational data 1 5 1 for galaxies with A" < 20 is shown as grey points. The upper 
panel shows the predictions from Bower et al. 1 3 1 and the lower one from the Baugh et al. 1 1 1 model. 
Histograms are normalised to give unit area under the distribution. Right: The median redshift as 
a function of K band magnitude for all galaxies (black) and for EROs with {R — K) colour redder 
than 5 (green), 5.3 (blue), 6 (red) and 7 (pink). For clarity, quartiles are only shown for EROs with 
[R — K) > 5.3. The upper panel shows predictions from Bower et al.|3| and the lower one from 
Baugh et al.| 1]. For comparison we plot the median redshift for all galaxies per magnitude bin in 
the A'-band observed in different studies (squares|2T|, diamonds |11 1 and circles |8|). The inset 
shows the median redshift per magnitude bin in the S-band for all galaxies predicted by Bower et 
al.|3 1 (solid line) and Baugh et al.| 1] (dashed line), and observational |9| as circles. 



4 Conclusions 

In this work we have extended the tests of the GALFORM galaxy formation code to 
include galaxies at z > 1 . We have tested two published models of galaxy formation, 
those of Baugh et al.yj and Bower et al.[3 1, against observed EROs number counts 
and K-selected redshift distributions. 

The Bower et al.[31 model gives an impressively close match to the number 
counts of EROs, while Baugh et al. 1 1 1 predicts an order of magnitude too few EROs. 
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Also, Baugh et al.Jl] model predicts a too low median redshift for K-selected galax- 
ies. We have experimented with Baugh et al.(|T| model without finding a way to 
reconcile its predictions with observations. The key seems to be that Bower et al.f3l 
model gives a better match to the observed evolution of the K-band luminosity func- 
tion, which means that this model puts massive galaxies in place earlier than in the 
Baugh et al.fl] model. This difference between the two models arises from the dif- 
ferent redshift dependence of the feedback processes which suppress the formation 
of massive galaxies and from the choice of the star formation time-scale. 
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